Non-Fermi liquid behavior and Griffiths phase in /-electron compounds 
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We study the interplay among disorder, RKKY and Kondo interactions in /-electron alloys. We 
argue that the non-Fermi liquid behavior observed in these systems is due to the existence of a 
Griffiths phase close to a quantum critical point. The existence of this phase provides a unified 
picture of a large class of materials. We also propose new experiments that can test these ideas. 
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The observation of non-Fermi-liquid (NFL) exponents 
in the thermodynamic and transport properties of /- 
electron alloys has stimulated considerable interest in the 
study of these materials |l| . The alloys in which NFL be- 
havior is observed fall into two categories: 1) Kondo hole 
systems, in which the /-electron atoms (R) are replaced 
by non-magnetic metallic atoms (M) according to the 
formula Ri^Mj, and 2) disordered ligand systems, in 
which the metallic atoms are substituted for a different 
metallic atom according to the formula R(M1)i_j,(M2)j,. 
Notice that due to alloying these compounds have a high 
probability of being disordered. That disorder is indeed 
a very important factor in bringing about the NFL be- 
havior in these compounds has been shown in recent ex- 
periments gj. This is in addition to the fact that most 
of these systems are close to a phase transition. Then, 
we claim that the NFL properties of these compounds 
are a consequence of the competition between the intra- 
site Kondo and the inter-site Ruderman-Kittel-Kasuya- 
Yosida (RKKY) interactions taking place in the midst of 
a disordered environment. If disorder were not present, 
there are two possibilities: the compound will have long- 
range magnetic order when the RKKY interaction is suf- 
ficiently large compared with the Kondo interaction, or 
the compound will be paramagnetic due to the quenching 
of the magnetic moments of the rare earth atoms. How- 
ever, the experimental observations show that the NFL 
behavior generally appears between these two phases [Q . 
Several proposals have been put forward as the possible 
explanations for the NFL behavior. A possible scenario 
is based on single impurity models with particular sym- 
metries such as multichannel Kondo effect of magnetic 
H and electric origin Another possible scenario 

attributes the NFL behavior to proximity to a quantum 
critical point Recently, a route to NFL behavior 

that emphasizes a disorder-driven mechanism, known as 
"Kondo disorder" , has been suggested . All of these 
proposals have had partial success in explaining some of 
the experiments In particular, conformal invariance 
scaling gives a good description of the dynamic suscep- 
tibility ~x"(lv, T) and the electrical resistivity in UCus-a; 
Pd x po|| ; a Kondo disorder model can explain the tem- 
perature dependence of the electrical resistivity ||; and 
the phase diagram of spin glass models is qualitatively 
similar to the one observed in these alloys M . Neverthe- 



less, there is no final word as to the origin of the NFL 
behavior in these alloys. Here, we propose a framework 
that incorporates what we believe are the essential as- 
pects of the problem: disorder and the competition be- 
tween RKKY and Kondo effect. In this framework the 
presence of disorder leads to the coexistence of a metallic 
(paramagnetic) phase with a granular magnetic phase. 
We show that this coexistence phase is equivalent to the 
Griffiths phase of dilute magnetic systems In our 

scenario, we have two electronic fluids: one of them is 
quenched by the Kondo interaction, behaving as a Fermi 
liquid; and the other is dominated by the RKKY interac- 
tion, leading to ordered regions of two level systems. We 
have, therefore, an inhomogeneous environment which is 
brought about by disorder. This scenario is reminiscent 
of the one found in compensated doped semiconductors 
(Si:P, Ge:Sb) p^-|l^]. In these systems, disorder leads to 
local density fluctuations which result in the formation 
of magnetic moments. The Griffiths phase is character- 
ized by the formation of rare strongly-coupled magnetic 
clusters which have large susceptibilities. In this phase 
the thermodynamic functions show essential singularities 
with strong effects at low temperatures. At these low 
temperatures, clusters of interacting magnetic moments 
can be thought of as "giant" spins which can tunnel over 
classically forbidden regions. In the Griffiths phase mag- 
netic clusters with N spins have a relaxation time which 
is given by [^|| (we use units such that K = ks = 1) 

T R = , (1) 

where luq is an attempt frequency and £ is a character- 
istic parameter that we discuss below. Due to cluster 
formation in the paramagnetic phase of these systems, 
we have the following predictions for the thermodynamic 
functions: 

7 = C v /T <x [ X (T)] av ex T- 1+x , 
[Xnl(T)] av c<T- 3 +\ 
[Xi(«0"L «^ 1+A tanh( W /T) , 

Tf^w) oc uj- 2+x Tta,nh(uj/T) , 
6 X (T)/x(T) ex T- A /2 , (2 ) 

where [...] av means average over disorder. Cy is the spe- 
cific heat; x(^) is the static susceptibility; Xni(T) is the 
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non-linear static susceptibility; xl(u) is the local fre- 
quency dependent susceptibility; 1/Ti is the NMR relax- 
ation rate; and <5x(T) is the mean square deviation of the 
susceptibility due to the distribution of susceptibilities in 
the system Jig . Here A = — £ _1 rn(c), where c denotes 
the density of the spins. The Griffiths phase is char- 
acterized by A < 1 so that the susceptibilities diverge at 
zero temperature. We propose that Griffiths singularities 
dominate the physics of the system at low temperatures 
leading to NFL behavior. Let us note here that power- 
law behaviors for 7 and \ have also been obtained by 
other researchers using different approaches 
Notice that the logarithmic behavior observed in some 
NFL compounds M can as well be fitted by small power 
laws (A w 1). Also, it follows from (|J) that NFL systems 
should have positively divergent non-linear susceptibili- 
ties (A < 3). Indeed, Uo.9Tho.iBei3 shows a tendency 
to a positively divergent susceptibility Et]], in contrast 
to the usual negative divergence of the paramagnetic 
susceptibility of UBei3. Systems like UCus_ a Pd y show 
even stronger divergent behavior and can be considered 
"deep" inside of the Griffiths phase. Recent neutron scat- 
tering experiments show that the imaginary part of the 
susceptibility, the specific heat and the static susceptibil- 
ity can be exactly fitted by the result (Q) with A = 2/3 
To compare this result with NMR and /iSR, we cal- 



culated the variation of the Knight shift, SK/ K cx 5\l 



for the same material 
for A = 2/3. 



Our result is shown in Fig 
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FIG. 1. Mean square deviation of the Knight shift as a 
function of temperature given in (Q). Diamonds: experiments 
for UCu 4 Pd from ref.[7]. 

The agreement between theory and experiment is very 
good. Our predictions are robust in the sense that all 
the results in (Q) have to be self-consistent. Moreover, 
we have a definite prediction for the NMR relaxation 
rate 1/Ti (a;), which should be largely frequency depen- 
dent. We also predict that under pressure the exponent A 
should change inside of the Griffiths phase Jil| . It would 
be interesting to plot the logarithmic derivative of the 
susceptibility (and/or specific heat) as a function of tem- 
perature and pressure, X(T,P) = 1 + dlnx(T, P)/d\nT, 
for various NFL systems to verify our predictions. The 



main characteristics of /-electron systems studied here 
are: interplay between RKKY and Kondo effects; mag- 
netic anisotropy; and disorder due to alloying. It is well 
known that /-electrons systems are characterized by their 
strong magnetic properties. This magnetism arises from 
the crystal field interactions and the strong spin-orbit 
coupling; in some of these materials, the large /- elec- 
tronic clouds lead to an anisotropy comparable in mag- 
nitude to the exchange energy. In addition to this ion 
anisotropy, one expects a Dzyaloshinsky-Moriya (DM) 
exchange interaction [ pp[ generated by the spin-orbit cou- 
pling. This type of interaction is only forbidden in highly 
symmetric situations which are not commonly realized in 
alloys with very complex unit cells . That this is the 
case, it can be seen from the coexistence of weak parasitic 
ferromagnetism within the antiferromagnetic phase in 
some materials such as R-Cr03 systems |22j| . We believe 
that this DM interaction accounts for the recent neutron 
scattering data in the NFL compound CeCu6-j,Auj, [ p3| . 
Consequently the RKKY and Kondo interactions will be 
strongly anisotropic in the /-electron systems. The sim- 
plest Hamiltonian that describes the situation above is 
the anisotropic Kondo model 



H = H e + Y^ 



JZ qZ Z 
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where the sum is over all the R atoms with spin S n . 
H e is the conduction electron Hamiltonian that can be 
obtained from band structure calculations, and s2 = 



E 



, r' 



with a — x,y,z is the conduction 



electron spin (r a are Pauli matrices). The main differ- 
ence between Kondo hole and disordered ligand systems 
is the alloying process. In the Kondo hole the magnetic 
atom is replaced which leads to a reduction of the number 
of magnetic moments. Moreover, the substitute atom has 
a different size from the original R atom which will lead 
to local changes in the lattice structure. In disordered 
ligand systems the lattice spacing is also affected by the 
substitution of the ligand atoms. As is well known, hy- 
bridization matrix elements between the conduction band 
and the /-electron system are exponentially sensitive to 
changes such as the type of lattice and characteristics of 
the substitute atom, and this will in turn affect the local 
values of the exchange constants in (||). Therefore, alloy- 
ing leads to a situation where local variations result in 
parts of the sample having larger exchange interactions 
than others. The magnetism in these systems can be un- 
derstood by looking at the effective interaction between 
the R atoms due to the conduction electrons. The mag- 
netic Hamiltonian in second order perturbation theory 
relative to the free electron problem is 



H 



M 
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(4) 



'' is the RKKY interaction mediated by the 
conduction electrons. In the limit of large anisotropy, 



where 
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J z » J- 1 , we have T z ' z ~ 0[(J z ) 2 /E F ] where E F is 
the Fermi energy. In a clean d-dimensional system this 
interaction decays like l/r d . In the presence of disorder, 
however, it decays like e~ r / e where t is the spin-orbit 
diffusion length |24| . Since we are dealing with disordered 
systems the effective interaction V% m is actually short- 
ranged. The important point here is that the ordering 
temperature of the magnetic system, T c , is of order of the 
exchange interaction, that is, T C (J Z ) cx T cx (J z ) 2 /Ep. 
The critical ordering temperature is not the only energy 
scale in the problem since the Kondo effect can take place 
and quench the magnetic moment. 

Consider for instance a particular position in the sys- 
tem, say n = M, where, due to the alloying, J| f is 
much larger than average. For simplicity we disregard 
all the other sites and look at the physics at this partic- 
ular site. The Hamiltonian of interest can be written as 
Hm = Hq + Hi, where 

H o = H e + Jm^u ( n M,-\ — n M,[) 



Sm s m) 



(5) 



and in the limit of large anisotropy we treat Hi as a 
perturbation of Hq. Observe that Hq can be easily diag- 
onalized in the S z basis. If the R atom is in a state such 
that S z is f|" (JJ.) the energy of the system is minimized by 
making a bound state with an electron with spin j (f), 
respectively. The bound state energy is just the Kondo 
temperature 



T K {M) = W e 



-1/(AT(0)J| 



(G) 



where iV(0) is the renormalized density of states at the 
Fermi energy, and W is the conduction electron band- 
width. This situation is similar to the approach to NFL 
behavior which takes into account a disordered distribu- 
tion of Kondo temperatures The ground state of 
Hq is doubly degenerate corresponding to the spin con- 
figurations | ff, |) and | t)- Application of first order 
perturbation theory in Hi shows that the singlet state 
I 1T> I) — I -II") T) is lower in energy than the triplet state 
I 1T> I) + I T) by an amount J^. In this case Hi acts as 
a transverse field and lifts the degeneracy of the bound 
state. In fact, using the bosonized version |2^] of Hamil- 
tonian © and mapping the problem into the dissipative 
two-level system pq] , we can show that tracing out the 
electrons of the problem results in a Sfj operator for the 
R atoms spin degrees of freedom [M . The competition 
between the RKKY interaction and the Kondo effect can 
be understood in terms of the two relevant energy scales, 
that is, T K {M, J z ) and T C {J Z ): i) if T K {M) » T c as we 
lower the temperature of the system below T K (M) the 
local magnetic moment is quenched and order is inhib- 
ited; ii) if Tk{M) «T c and the temperature is lowered 
below T c there is local magnetic order and the Kondo 
effect is suppressed [^8] . If we now take into account the 
magnetic moments that are left to interact via RKKY, 



as in (||), together with the sites which are quenched by 
the Kondo effect, as in (|J), we see that the magnetism 
of the original problem in the limit where J„ >> is 
described by 



H, 



eff 



Ti t jSl Sj + tiS i 



(7) 



where the brackets imply nearest neighbor coupling. 
Notice that Tij ~ 0{{J z ) 2 /E F ) and U ~ ©(J- 1 ) are ran- 
dom (but inter-correlated) variables dependent on the 
alloying and lattice structure. Thus we have mapped 
our problem into the random Ising model in a random 
transverse magnetic field [[l5| . The phase diagram of this 
model follows: at small doping the RKKY interaction 
dominates and the system can order magnetically (the 
ordered phase can be antiferromagnetic, ferromagnetic or 
spin glass j7| depending if the mean value of r,-^ in (^) is 
positive, negative or null, respectively). With increasing 
doping the quantum fluctuations grow due to the Kondo 
effect and the bulk critical temperature decreases until it 
vanishes for some critical value of doping. At this quan- 
tum critical point the system percolates. For large values 
of doping, inside of the paramagnetic phase, only finite 
clusters of magnetic atoms can be found. Among these 
clusters there are some rare ones which are large and 
strongly coupled. Within these clusters the spins behave 
coherently as a "giant" spin or a magnetic grain. Wc 
can describe the cluster in terms of a "classical" degree 
of freedom which can be parametrized by Euler angles 
(0,(p). The classical energy, E(9,<f>), will have at least 
two minima due to the original degeneracy of the mag- 
netic ground state. In the simplest of cases, E(9, <fi) can 
be written in terms of a classical spin with X-easy axis 
and XF-easy plane, 



E(9, ( f>) = N(-e x + e ll 



sin 



sin 



(8) 



where e± > en > are the anisotropy energies per- 
pendicular and parallel to the easy axis, respectively. 
These energies depend on the microscopic coupling con- 
stants in (mj. Observe that the energy has two minima 
at (tt/2,0) and (tt/2,tt) with an energy barrier between 
them. When the temperature is higher than the bar- 
rier height the cluster is thermally activated and behaves 
classically. At lower temperatures the cluster can un- 
dergo quantum tunnelling between the two minima. Us- 
ing standard instantons methods |30| , we can calculate 
the parameters that appear in (m): 



uj = 2^ejJeT 



C w In 



1 + aA||A± 



(9) 



Notice that in the absence of anisotropy (e|| = ej_) tun- 
nelling cannot occur (£ — > oo and — > oo), as expected. 
Then, the low energy physics for a cluster £1 with N spins 
reduces to the Hamiltonian 
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H 



where A^at) = 1/tr is the tunnelling energy and is given 
in M), and it can be related to the anisotropy energies 
by (§). Using @) and averaging over cluster with dif- 
ferent sizes, we arrive at the predictions given in (|^) 
p6[ . In conclusion, we propose that the NFL behav- 
ior observed in /-electron systems can be attributed to 
the existence of Griffiths singularities close to the quan- 
tum critical point. These singularities have their origin 
on the interplay between the RKKY and Kondo interac- 
tions in the presence of magnetic anisotropy and disorder. 
These conclusions have similarities with those discussed 
recently by Sachdev (^9|. We were also able to map the 
disordered Kondo lattice problem into the random Ising 
model in a random transverse magnetic field where dis- 
order is correlated. In the paramagnetic phase of this 
Ising model the physics of clusters can be understood 
in terms of the quantum tunnelling of intrinsic magnetic 
grains which are described by a classical spin model. At 
low temperatures the spin degrees of freedom of a mag- 
netic grain can tunnel over classically forbidden regions 
and at finite temperatures they can be thermally acti- 
vated. At very low temperatures the problem reduces to 
a two-level system problem which, when appropriately 
averaged over disorder, leads to the Griffiths singulari- 
ties and to the predictions in (||). This Griffiths phase 
will depend strongly on the type of lattice structure and 
value of the local microscopic exchange constants. This 
would explain why systems like CeCu6- y Au y have to be 
fine tuned for NFL behavior to be observed, while others, 
like Thi_ x U x Pd3, have large regions of NFL behavior. It 
is indeed possible, as in ID systems [p|, that Griffiths 
behavior can extend over large regions of doping |jT| . 
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Note added.- After this paper was completed, exper- 
imental indications of possible Griffiths phase behavior 
have been reported 133] . 



[1] See, for instance, Conference on Non-Fermi liquid behav- 
ior in Metals, J. Phys. Cond. Mat 8, 1996. 

[2] T. Graf et al, Phys. Rev. Lett. 78, 3769 (1997). 

[3] P. Nozieres and A. Bladin, J. Physique 41, 193 (1980); 
B. Andraka and A. M. Tsvelik, Phys. Rev. Lett. 67, 2886 
(1991); I. Affleck and A. W. W. Ludwig, Nucl. Phys. 
352, 849 (1991); P. D. Sacramento and P. Schlottmann, 



(10) Phys. Rev. B 43, 13294 (1991). 

[4] D. L. Cox, Phys. Rev. Lett. 59, 1240 (1987). 
[5] Studies of the lattice effects in these models have just 
begun: F. Anders et al, Phys. Rev. Lett. 78, 2000 (1997). 
[6] J. A. Hertz, Phys. Rev. B 14, 1165 (1976); A. J. Mil- 
lis, ibid. 48, 7183 (1993); A. M. Tsvelik and M. Reizer, 
ibid. 48, 9887 (1993); S. Sachdev and J. W. Ye, 
Phys. Rev. Lett. 69, 2411 (1992). 
[7] A. M. Sengupta and A. Georges, Phys. Rev. B 52, 10295 

(1995); S. Sachdev and N. Read in Ref. 
[8] O. O. Bernal, et al., Phys. Rev. Lett. 75, 2023 (1995). 
[9] V. Dobrosavljevkf et al., Phys. Rev. Lett. 69, 1113 
(1992); E. Miranda, et al, ibid. 78, 290 (1997). 
[10] M.C. Aronson, et al., Europhys. Lett., 40, 245 (1997). 
[11] R. B. Griffiths, Phys. Rev. Lett. 23, 17 (1969). 
[12] R. N. Bhatt and P. A. Lee, Phys. Rev. Lett. 48, 344 

(1982); M. Milovanovic et al, ibid. 63, 82 (1989). 
[13] A. Langenfeld and P. Woffle, A. Phys. 4, 43 (1995). 
[14] R. N. Bhatt and D. Fisher, Phys. Rev. Lett. 68, 4141 
(1992). This work was among the first to predict that 
disorder gives divergent Xi an d could lead to NFL be- 
havior in heavy fermions. 
[15] M. Guo, et al, Phys. Rev. Lett. 72, 4137 (1994); 
H. Rieger and A. P. Young, ibid. 72, 4141 (1994); 
D. Fisher, Phys. Rev. B 51, 6411 (1995); M. Guo et al., 
Phys. Rev. B 54, 3336 (1996). 
[16] For a derivation of some of these results, see for instance, 
M. J. Thill and D. A. Huse, Physica A 214, 321 (1995). 
[17] F. G. Aliev, et al, Physica B 206, 454 (1995). 
[18] M. C. Aronson, et al, Phys. Rev. Lett. 75, 725 (1995). 
[19] B. Bogenberger and H. v. Loneysen, Phys. Rev. Lett. 
74, 1016 (1995); K. Umeo,et al., J. Phys. Cond. Mat. 
48, 9743 (1996). 
[20] I. Dzyaloshinsky, J. Phys. Chem. Solids 4, 241 (1958); 

T. Moriya, Phys. Rev. 120, 91 (1960). 
[21] K. Yosida, Theory of Magnetism (Springer, Berlin, 1991). 
[22] W. C. Koeler, et al, Phys. Rev. 118, 58 (1960). 
[23] A. Rosch, et al., Phys. Rev. Lett. 79, 159 (1997). 
[24] A. Jagannathan, et al., Phys. Rev. B 37,436 (1988). 
[25] V. J. Emery and S. A. Kivelson, in Fundamental Prob- 
lems in Statistical Mechanics, VIII, H. van Beijeren and 
M. H. Ernst, eds. (Elsevier, Amsterdam, 1994), pg. 1; G. 
Kotliar and Q. Si, Phys. Rev.B 53, 12373. 
[26] A. J. Leggett et al., Rev. Mod. Phys. 59, 1 (1987). 
[27] This picture is valid for JljN(0) > 0.6, which has been 
reported (F.Lesage and H. Saleur, Phys. Rev. Lett. 80, 
4370 (1998)). 

[28] A model that explores this competition is the two impu- 
rity Kondo model:B. A. Jones, et al., Phys. Rev. Lett. 
61, 125 (1988). See also S. Doniach, Physica B & C, 91, 
231 (1977). 

[29] S. Sachdev, Phil. Trans. R. Soc, 356, 173 (1998). 
[30] N. V. Prokofev and P. C. E. Stamp, J. Low Temp. Phys. 

104, 143 (1996). 
[31] In Thi_a;Ua;Ru2Si2, NFL behavior is seen below x ~ 
0.07. See, H. Amitsuka et al, Physica B 186 337 (1993). 
[32] M. C. de Andrade, et al, |cond-mat/9802081. 



4 



